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L  Statement  of  Work 


Generally  a  missile  is  called  a  cruise  missile  if  its  speed  is  subsonic,  if  it  uses  a  built-in 
global  positioning  navigation  system  (GPS/INS),  and  if  its  range  is  at  least  several 
hundred  miles.  A  specific  mission  for  a  cruise  missile  is  programmed  by  specifying  a 
sequence  of  waypoints  and  including  any  TERCOM  maps  that  are  available.  The 
missile  uses  its  on-board  computational  facilities  and  its  GPS/INS  system  to  guide  it  to 
its  target. 

During  the  1990s,  cruise  missiles  launched  from  ships  and  submarines  were  used  with 
great  success  in  the  Persian  Gulf.  Targets  which  may  be  hundreds  of  miles  inland  have 
been  successfully  attacked  by  these  Navy  weapons  stationed  in  relatively  safe  positions 
offshore.  In  addition,  cruise  missiles  have  been  used  to  disable  SAM  sites,  making 
manned  aircraft  strikes  safer. 

By  representing  threat  regions  as  circles,  the  problem  of  developing  a  single  mission  for 
a  single  cruise  missile  can  be  viewed  as  a  computational  geometry  problem  in  which  we 
seek  a  path  composed  of  line  segments  from  a  launch  site  to  a  target  site  which  skirts 
selected  threat  regions.  Using  current  technology,  missions  are  developed  using  a 
two-step  process.  In  the  first  step,  a  mission  planner  uses  a  two-dimensional  map  to 
manually  select  a  path  from  a  launch  site  to  the  target.  In  the  second  step  a  software 
system  includes  the  vertical  dimension  and  develops  an  estimate  of  the  probability  of  a 
successful  mission.  This  process  may  be  repeated  several  times  until  the  mission 
planner  is  satisfied  with  the  mission  plan.  In  this  research,  we  present  three  algorithms 
designed  to  help  automate  the  first  step  of  the  mission  planning  process. 


11.  Finding  Safe  Paths 


In  this  study,  we  lay  a  grid  over  the  military  theatre  of  interest  and  create  a  grid  graph 
[N,  A].  The  points  s  e  N  and  t  e  N  are  defined  to  be  the  grid  points  nearest  the  missile 
and  target,  respectively.  For  each  edge  we  define  a  traversal  probability.  Traversal 
probabilities  within  a  threat  region  are  set  to  values  which  are  proportional  to  the 
distance  from  the  midpoint  of  an  edge  to  the  center  of  the  threat.  Edges  which  are 
outside  the  threat  regions  have  traversal  probabilities  very  close  to  1.  By  applying  a 
slight  modification  of  Dijkstra’s  algorithm,  we  can  easily  find  the  safest  path  from  s  to  t. 
The  set  of  nodes  which  appear  in  near— optimal  safe  paths  has  been  called  the  safe 
corridor.  This  manuscript  (which  appears  in  Appendix  A)  presents  a  proof  of 
correctness  for  both  the  problem  of  finding  the  safest  path  and  the  problem  of  finding  a 
safe  corridor  in  a  grid  graph. 


III.  Grid— Free  Algorithms 


In  this  study,  we  present  a  new  simple,  yet  surprisingly  effective  algorithm  that 
automatically  finds  a  path  from  the  launch  site  to  the  target  site  which  avoids  all  threats 
and  does  not  require  the  use  of  a  grid.  The  circumscribed  triangle  algorithm  constructs 
a  mission  by  moving  from  a  given  waypoint  tangentially  to  the  circumference  of  nearby 
threat  regions.  These  directions  are  determined  by  constructing  a  circumscribed 
triangle  around  a  threat  region.  This  basic  strategy  is  embedded  within  an  efficient 
branch -and -bound  framework.  On  twenty  randomly  generated  test  problems,  this 
algorithm  worked  extremely  well  in  all  cases.  This  study  is  documented  in  Appendix  B. 


rV.  Mission  Planning  With  A  Probability  Side  Constraint 


After  consultation  with  analysts  at  the  Naval  Surface  Warfare  Center  at  Dahlgren, 
Virginia,  we  discovered  that  real  mission  plans  frequently  penetrate  enemy  threat 
regions.  In  this  investigation,  we  extended  the  circumscribed  triangle  algorithm  to 
incorporate  a  side  constraint  on  the  probability  of  a  successful  traversal.  This  allows  the 
missile  to  take  short  cuts  through  threat  regions,  as  long  as  the  probability  of  success 
meets  some  user— specified  criteria.  This  work  is  documented  in  Appendix  C. 
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Abstract 

The  problem  of  constructing  strike  plans  for  cruise  missiles  motivated  the 
formulation  of  two  new  graph  problems:  the  safest  path  problem  and  the  safe 
corridor  problem.  The  idea  is  to  construct  a  flight  path  for  a  cruise  missile 
through  a  military  theatre  which  results  in  a  high  probability  of  mission 
success.  Exact  algorithms  for  both  problems  are  presented  along  with  a  proof 
of  correctness.  The  notion  of  a  safe  corridor  is  illustrated  on  several  randomly 
generated  strike  planning  problems  involving  a  single  cruise  missile  and  a 
single  target  protected  by  several  surface-to-air  missile  sites.  The  network  is 
deflned  by  a  grid  graph  over  the  terrain  of  interest. 
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1  Introduction 


Let  G  —  [iV,  A]  be  a  directed  graph  consisting  of  a  set  N  of  nodes  and  a  set 
A  of  arcs  whose  elements  are  ordered  pairs  of  distinct  nodes.  A  directed  path 
in  G  from  node  s  to  node  t  is  an  alternating  sequence  of  nodes  and  arcs  from 
G  of  the  form  P,t  =  {s  =  =  t}. 

Let  pij  denote  the  probability  of  a  safe  traversal  over  the  arc  (i,j).  Then 
the  probability  of  a  safe  traversal  over  the  directed  path  Pst  is  given  by 
PhhPi2i3---Pik-iik-  Given  s  €  iV  and  t  E  N,  the  safest  path  problem  is  to  find 
a  directed  path  in  G  having  the  largest  probability  of  safe  traversal.  Given 
R  <  I,  s  E  N  and  t  E  N,  the  safe  corridor  problem  is  to  find  the  set  of  nodes 
C  Q  N  such  that  for  each  n  E  C  there  exists  a  path  containing  n  whose 
traversal  probability  is  at  least  i?p*,  where  p'  is  the  traversal  probability  of 
the  safest  path  in  G. 

The  objective  of  this  investigation  is  to  develop  an  eflicient  algorithm  for 
the  safest  corridor  problem  and  demonstrate  how  this  can  be  used  to  plan 
flight  paths  for  cruise  missiles.  The  safe  corridors  generated  allow  a  mission 
planner  to  develop  flight  paths  which  maximize  mission  success. 

The  safest  path  problem  is  closely  related  to  the  classical  shortest  path 


problem  which  is  discussed  in  numerous  books  (see  for  example  [1,  2,  3,  7, 
10,  11, 12, 13,  14, 18]).  Excellent  surveys  of  the  many  shortest  path  problem 
variations  may  be  found  in  [5,  8,  9).  The  idea  of  a  safe  corridor  was  first 
introduced  to  the  authors  by  Solka  and  Rogers  [15],  who  were  developing  a 
safe  corridor  using  Dijkstra’s  shortest  path  algorithm.  Parallel  versions  of 
their  algorithms  may  be  found  in  [16,  17]. 

2  The  Safest  Path  Algorithm 

The  safest  path  algorithm  presented  in  this  section  is  a  straight-forward  adap¬ 
tation  of  Dijkstra’s  classical  algorithm  [6]  for  the  one-to-one  shortest  path 
problem.  The  input  is  a  directed  graph  [N,  A]  with  node  set  iV  =  {1, 2, ...,  n} 
and  arc  set  A.  Associated  with  each  arc  (i,j)  is  the  probability  of  a  suc¬ 
cessful  traversal,  p,j.  A  safest  path  is  desired  from  the  origin  s  €  AT  to  the 
destination  t  €  iV,  and  it  is  assumed  that  a  path  from  s  to  t  exists. 

2.1  The  Algorithm 

The  safest  path  algorithm  begins  at  node  s  and  constructs  a  safest  path  tree 
T  in  which  the  safest  path  from  s  to  any  node  in  T  is  known.  When  node  t 


is  placed  in  the  tree  we  have  a  safest  path  from  s  to  t  and  the  algorithm  ter¬ 
minates.  The  algorithm  uses  two  working  arrays  and  working  sets  Sk  which 
contain  the  nodes  in  the  current  T  at  iteration  k.  The  array  dj  denotes  the 
probability  of  traversal  from  s  to  j  and  the  array  bj  denotes  the  predecessor 
of  node  i.e.  if  (i,j)  6  T,  then  bj  -  i.  At  each  iteration,  T  is  enlarged  by 
one  node,  i.e.  \Sk\  =  k.  Hence,  in  at  most  \N\  iterations  the  algorithm  will 
terminate. 

Algorithm  :  safest  path 

Inputs  ;  N,A,s,t,pij 

Outputs  :  r,  bj  for  all  nodes  j  €  T,  and  dt 

Assumption  :  There  exists  a  path  from  s  to  t  in  [iV,  A]. 

begin 

Sq  < —  0;  k  < —  Oj 

dj  * - 1  for  all  j  6  iV; 

dg  < —  1|  bg  i —  0| 

while  t  ^  Sk  do 
begin 


select  ijt  €  N\Sk  such  that  Vj  €  N\Sk,  d,*  >  dj] 


Sk+1  <—  5fc  U  {ijt};  k  *—  k  +  1] 

for  all  j  €  N\Sk  such  that  (ijtji)  €  A  do 


begin 

if  dj  <  ^ikPiki  then  dj  <  ^ikP^kJi  *  *^5 
end; 
end; 
end; 

An  algorithm  for  the  one-to-all  safest  path  problem  can  be  obtained  by  mod¬ 
ifying  the  termination  criterion  to  “while  N\Sk  r  0 


2.2  Correctness  of  the  Algorithm 

^  In  this  section  we  show  that  the  safest  path  algorithm  terminates  with  the 

correct  solution.  The  proof  is  based  on  induction  and  relies  on  the  fact  that 
the  safest  path  from  s  to  every  node  in  the  safest  path  tree  is  known. 

• 

Theorem.  When  the  safest  path  algorithm  terminates,  dt  is  the  probability 
of  the  safest  path  from  s  to  f  in  [N,  A]. 

Proof.  We  proceed  by  finite  induction  using  the  predicate 
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P{m):  (a)  for  all  i  E  Sm,  di  is  optimal  with  an  optimal  path 


from  s  to  i  lying  in  Sm,  and 
(b)  for  all  j  €  N\Sm  with  dj  >0,  dj  is  the  largest 

traversal  probability  for  a  path  from  s  to  j  in  which 
all  nodes  except  j  lie  in  Sm,  and  such  a  path  can  be 
decomposed  into  a  subpath  from  s  to  bj  and  the 
subpath  consisting  of  bj  followed  immediately  by  j. 

When  5o  =  0,  only  dg  >  0.  So  Si  =  {s}  and  dj  =  1.  Since  any  path  with 
probability  1  must  be  optimal,  dg  is  optimal  and  the  degenerate  optimal  path 
from  s  to  s  lies  in  5i.  Thus  (a)  holds.  Also,  the  only  path  from  s  to  a  node 
j  in  N\Si  consists  of  s  and  the  arc  connecting  s  to  jf.  For  any  such  point  j, 
dj  was  set  to  dgpgj  =  pgj,  the  traversal  probability  associated  with  the  only 
path  from  the  single  node  s  of  to  j,  using  only  bj  =  s  and  j.  So  (b)  holds. 
Thus  P(l)  is  true. 

Let  k  be  an  arbitrary  fixed  integer  such  that  1  <  k  <  n.  Assume  that 
P(l), . . . ,  P{k)  are  all  true.  Consider  Sk+i  and  ik  €  5fc+i\5fc.  Since  P{k)  is 
true,  di  is  optimal  for  all  i  E  Sk,  and  optimal  paths  from  s  to  those  nodes 
also  lie  in  Sk-  Let  Pq  be  the  path  from  s  to  ik,  with  traversal  probability  po, 


consisting  of  the  path  from  s  to  6,-^  followed  immediately  by  ik- 

Suppose  (a)  does  not  hold  for  Sk+i-  It  must  be  that  there  is  a  path  Pi 
from  s  to  ik  with  traversal  probability  pi  >  dij. .  Furthermore,  this  path  must 
contain  at  least  one  node  of  N\Sk+i  •  Let  q  be  the  first  node  of  N\Sk+i  on 
path  Pi,  so  that  Pi  decomposes  into  two  subpaths  P2  from  s  to  5  and  P3 
from  q  to  ik.  Let  p2  eind  pz  be  the  traversal  probabilities  of  paths  P2  and  P3, 
respectively.  All  nodes  of  P2  except  q  lie  in  Sk-  So  by  (b),  d,  is  the  largest 
traversal  probability  for  a  path  from  s  to  g  in  which  all  nodes  except  q  lie  in 
Sk-  But  when  ik  was  selected,  d,-,.  >  d,.  So  pi  =  pzpz  <  d,p3  <  d,,.p3  <  d,-,^, 
a  contradiction.  Thus  (a)  holds  for  Sk+i- 

We  now  show  (b)  holds  for  Sk+i-  Let  r  €  N\Sk+i  with  dr  >  0.  Let  P4 
be  the  path  with  the  largest  traversal  probability  for  a  path  from  s  to  r  in 
which  all  nodes  except  r  lie  in  Sk+i  and  let  that  probability  be  p4.  Let  u  be 
the  last  node  of  Sk+i  on  P4. 

Case  1.  Suppose  that  before  ik  was  selected,  dr  =  —1.  Then  after  ik  was 
selected,  the  only  paths  from  5  to  r  with  nodes  except  for  r  lying  entirely  in 
Sk+i  have  ik  and  r  as  the  last  two  nodes  and  (b)  holds  for  r. 

Case  2.  Suppose  that  before  ik  was  selected,  d^  >  0.  Thus  before  ik  was 
selected,  since  (b)  holds  for  Sk,  dr  was  the  largest  traversed  probability  for  a 
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path  from  5  to  r  in  which  all  nodes  except  r  lie  in  Sk-  Let  pr  be  the  value  of 
dr  at  that  point. 

Subcase  2.1.  Assume  u  is  4.  Then  =  di^pi^r-  Thus  following  the  ad¬ 
justment  to  dr  after  ijt  was  selected  dr  =  p4  since  p4  is  the  largest  probability 
path  and  (b)  holds  for  Sk+i- 

Subcase  2.2.  Assume  u  is  not  ik,  but  is  some  other  node  j  G  Sk^i-  Thus 
P4  will  not  contain  ik  and  no  adjustment  was  made  to  dr  after  4  was  selected, 
so  that  dr  =  Pr  and  (b)  holds  for  5fc+i. 

I 


3  The  Safe  Corridor  Algorithm 

For  applications  in  the  area  of  strike  planning,  there  are  many  alternate 
optima.  In  addition,  there  are  many  paths  that  come  very  close  to  being 
optimal.  The  set  of  nodes  which  appear  in  near-optimal  safe  paths  is  Cctlled 
the  safe  corridor. 

The  safe  corridor  algorithm  constructs  two  safest  path  trees,  Ti  and  T2- 
T\  records  .the  safest  path  from  s  to  any  node  in  its  tree,  and  T2  records 
the  safest  path  from  t  to  any  node  in  its  tree.  All  nodes  in  the  graph  are 


initially  assigned  a  negative  probability  of  successful  axrival  from  the  source 
and  from  the  destination.  Some  nodes  may  not  be  updated  with  non-negative 
probabilities,  but  all  nodes  near  a  safest  path  will  be  updated. 

For  any  node  j,  let  6j  denote  the  greatest  probability  of  successful  traver¬ 
sal  from  node  j  to  the  destination,  t,  let  Pj  denote  the  predecessor  of  node 
j  in  T2,  and  let  (k  denote  the  set  containing  the  current  nodes  in  T2  at  it¬ 
eration  k.  The  greatest  probability  of  successful  traversal  from  the  source 
to  the  destination  via  node  j  is  djSj,  where  dj  and  Sj  are  the  safest  route 
probabilities  from  s  to  j  and  from  j  to  t,  respectively,  and  have  been  updated 
with  nonnegative  probabilities.  For  any  node  i  along  the  safest  path  from 
s  to  t,  diSi  =  dt.  Also  Ss  =  dt,  since  pij  =  pj,-.  H  djdjjdt  is  very  close  to  1, 
then  node  j  is  not  on  a  safest  path,  but  the  probability  of  successful  traversal 
from  s  to  t  via  j  is  almost  as  great  as  that  of  a  safest  path.  Hence,  a  safe 
corridor  can  be  generated  by  including  all  nodes  j  for  which  dj6j  is  within  a 
given  range  near  dt. 

The  safe  corridor  C  C  N  is  determined  by  a  pairameter  i?  <  1.  K  i2  =  1, 
then  every  node  in  C  is  on  some  optimal  path  from  s  to  t.  If  i?  =  0.9,  then 
the  nodes  in  C  are  on  paths  whose  traversal  probability  is  at  most  10%  less 
than  the  probability  of  the  safest  path.  If  i?  =  0,  then  C  =  N.  The  algorithm 


is  given  as  follows: 


Algorithm  :  safe  corridor 
Inputs  :  N,A,s,t,pij,R 
Outputs  :  C 
begin 

Sq  < —  05  k  i —  O5 

dj  i - 1,  6j  i - 1  for  all  j  €  N] 

d$  1;  O5  5t  ^1;  A  ^  0; 

while  t  ^  Sk  do 

begin 

select  ifc  €  N\Sk  such  that  Vj  €  N\Ski  ^  dj", 
Sk+i  <—  5jt  U  {ifc};  fc  <—  +  1; 
for  all  j  €  N\Sk  such  that  (zfc,  j)  6  A  do 
begin 

if  dj  <  dij^pii^j  then  dj  <  d\^pif^j\  bj  <“  it, 
end; 
end; 

Co  <-  0;  k*-Q\ 
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while  3  ^  Cit 
begin 

select  hk  G  iV\Cfc  such  that  Vj  €  N\^k, 

Cjt+i  Cfc  U  {hk}]  k  *—  k  +  1] 

for  all  j  G  N\(k  such  that  (j,  hk)  G  A  do 

begin 

if  6j  <  hkPjhk  then  6j  <-  Sh^pjh),]  *-  hk] 
end; 
end; 

for  all  j  G  iV 
begin 

if  djSj  >  R{dt)  then  C  C  U  {;}; 
end; 
end; 

The  value  of  R  must  be  carefully  chosen  for  each  problem.  A  value  too 
close  to  1  can  exclude  valuable  options.  However,  as  the  value  of  R  is  reduced, 
the  safe  corridor  can  very  quickly  expand  to  include  a  large  portion  of  the 
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4  The  Two-Dimensional  Strike  Planning 
Problem  on  a  Grid  Graph 


Strike  planning  for  cruise  missiles  involves  finding  a  flight  path  through 
a  combat  theatre  that  involves  threats.  The  most  serious  threat  for  a  cruise 
missile  is  an  air  defense  system  consisting  of  a  search  radar,  a  fire-control 
radar,  and  a  SAM  weapon  (see  [19]).  The  effectiveness  of  the  SAM  system  - 
depends  on  the  total  amount  of  time  it  has  to  engage  and  destroy  the  cruise 
missile.  The  search  radar  is  usually  always  active,  while  the  fire-control 
radar  only  becomes  active  when  needed.  The  search  radar  spends  some  time 
in  determining  i|  a  missile  detect  is  real  or  is  a  false  alarm.  A  cruise  missile 
flies  at  low  altitude  in  an  attempt  to  avoid  detection.  The  fire-control  radar 
is  turned  on  only  after  search  radar  personnel  decide  that  they  have  a  real 
target.  There  is  some  warmup  required  and  some  time  needed  to  aim  and 
Are  the  SAM.  The  less  time  a  cruise  missile  spends  near  one  of  these  sites, 
the  higher  the  probability  of  mission  success. 


For  this  application,  we  lay  a  grid  over  the  area  and  form  a  grid  graph 
[iV,  A].  The  points  s  ^  N  and  f  €  iV  are  the  grid  points  nearest  the  missile 
and  target,  respectively.  The  traversal  probabilities  are  related  to  the  loca¬ 
tion  of  the  threats.  Given  R  <  1,  the  objective  is  to  determine  C  Q  N,  the 
safe  corridor  through  [iV,  A]. 

To  demonstrate  the  safe  corridor  algorithm,  we  developed  a  random  prob¬ 
lem  generator.  Each  interior  point  in  the  grid  has  degree  8  as  illustrated  in 
Figure  1.  The  missile  is  randomly  placed  on  the  left  vertical  boundary  and 
the  target  is  randomly  placed  on  the  right  vertical  boundary.  The  SAM 
systems  are  placed  at  random  nodes  in  [A*.  A].  Three  test  problems  and 
the  corresponding  safe  corridors  are  illustrated  in  Figures  2,  3,  and  4.  The 
traversal  probabilities  are  in  the  interval  [0.1.  0.99].  Probabilities  within  a 
threat  region  (determined  by  the  distance  to  the  fire-control  radar)  are  pro¬ 
portional  to  the  distance  from  the  midpoint  of  an  edge  to  the  fire-control 
unit.  Horizontal  and  vertical  edges  which  are  outside  threat  regions  have 
traversal  probabilities  of  0.99.  Diagonal  edges  have  traversal  probabilities  of 
0.985.  Hence,  traversing  a  diagonal  arc  is  slightly  safer  than  traversing  one 
horizontal  and  one  vertical  arc  connecting  the  endpoints  of  the  diagonal  arc. 


Figures  1  through  4  about  here 


5  Summary  and  Extensions 

This  manuscript  presents  a  new  graph  problem  which  we  call  the  safest  path 
problem.  A  Dijkstra-like  algorithm  for  this  problem  is  presented  along  with 
a  correctness  proof.  This  work  was  motivated  by  the  problem  of  developing 
strike  plans  and  we  demonstrate  how  our  algorithm  can  be  used  to  develop 

safe  corridors  on  two-dimensional  grid  graphs. 

The  work  presented  in  this  manuscript  can  be  extended  to  three  dimen¬ 
sions.  Even  though  the  algorithm  is  very  fast,  problem  size  could  eventually 
lead  to  excessive  computation  time.  It  is  also  the  case  that  repeated  flights 
over  a  threat  may  cause  threat  personnel  to  become  more  alert  so  that  the 
second  flight  near  a  threat  may  have  a  lower  probability  of  successful  traver¬ 
sal  than  the  first.  Planning  for  multiple  strikes  would  need  to  account  for 
this.  A  description  of  this  problem  along  with  solution  approaches  may  be 


found  in  [4,  20]. 
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Abstract 


The  problem  of  mission  planning  for  a  cruise  missile  in  two  dimensions  can 
be  viewed  as  a  mathematical  problem  in  the  area  of  computational  geometry 
in  which  we  seek  a  path  composed  of  a  set  of  line  segments  from  a  given 
missile  location  to  a  given  target  location  which  does  not  pass  through  any 
SAM  radar  site.  The  shortest  path  requiring  the  fewest  number  of  turns  is 
most  desirable.  Two  algorithms  for  this  problem  have  been  developed  and 
compared  in  an  empirical  study. 
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1  INTRODUCTION 


For  this  investigation,  we  clcissify  a  missile  as  a  cruise  missile  if  its  speed  is 
sub-sonic,  it  uses  a  built-in  global  positioning  navigation  system  (GPS/INS), 
and  its  range  is  at  least  300  miles.  A  specific  mission  for  a  cruise  missile 
is  programmed  by  simply  specifying  a  sequence  of  coordinates.  The  missile 
uses  its  on-board  computational  facilities  and  its  GPS/INS  system  to  guide  it 
through  this  sequence  of  points.  At  present,  once  a  cruise  missile  is  launched, 
its  mission  cannot  be  modified. 

A  cruise  missile  can  be  destroyed  in  flight  by  a  SAM  air  defense  system 
consisting  of  three  components:  a  search  radar,  a  fire-control  radar,  and 
a  SAM  weapom.  A  complete  description  of  this  system  may  be  formd  in 
Zuniga  and  Gorman  [7].  According  to  Zuniga  and  Gorman,  the  effectiveness 
of  the  SAM  system  depends  on  the  total  amount  of  time  it  has  to  engage 
and  destroy  the  cruise  missile.  The  search  radar  only  becomes  active  when 
needed.  Some  time  is  spent  to  determine  if  a  missile  detect  is  real  or  false. 
If  the  detect  is  determined  to  be  real,  then  the  fire-control  radar  is  turned 
on  to  aim  and  fire  the  SAM.  The  search  radar  has  the  same  effective  range 
which  is  represented  by  a  circle  in  two  dimensions. 


Mission  planning  for  a  single  cruise  missile  involves  constructing  a  flight 
path  from  a  given  origin  to  a  given  destination  through  a  combat  theatre 
having  SAM  sites.  In  two  dimensions,  the  SAM  sites  are  represented  as 
circles  and  the  objective  is  to  find  a  set  of  line  segments  linking  the  origin  to 
the  destination  which  do  not  pass  through  any  SAM  threat  circle. 

Most  research  investigations  in  the  area  of  mission  planning  begin  by  plac¬ 
ing  a  grid  over  the  combat  theatre  and  applying  a  modification  of  Dijkstra’s 
algorithm  (see  [2])  to  obtain  a  safe  path  from  the  origin  to  the  destination 
(see  [1,  3,  4,  5,  7,  8]).  There  are  two  disadvantages  to  this  approach.  The  grid 
graph  can  become  very  large,  particularly  for  the  three-diniensional  case,  and 
consequently  the  algorithms  can  be  very  time  consuming.  Also,  the  resulting 
paths  may  involve  too  many  line  segments.  A  path  with  only  a  few  line  seg¬ 
ments  is  viewed  as  better  than  one  with  many  line  segments.  The  objective 
of  this  investigation  is  to  present  new  algorithms  for  mission  planning  that 
do  not  require  the  use  of  a  grid. 


2  THE  MISSION  PLANNING  PROBLEM 


The  problem  of  finding  a  feasible  mission  plan  can  be  defined  mathematically 
as  follows; 

Given  two  points  M  and  T  G  and  K  circles  each  with  center 
Xi, Xfc  €  Ef  and  radii  ri, tk",  find  a  set  of  line 
segments  [¥,,¥2],  [¥s,¥s+,]  with  Tj  =  M  and 

Ts+i  =  T  such  that  no  line  segment  intersects  any  circle. 

A  feasible  mission  plan  through  a  theatre  having  five  SAM  sites  is  illustrated 
in  Figure  1.  The  quality  of  a  mission  plan  is  measured  by  the  length  of  the 
path  and  the  number  of  segments.  The  best  path  is  the  straight  line  from 
M  to  T.  If  this  is  not  feasible,  then  we  seek  a  short  path  with  only  a  few 
segments. 


Figure  1  about  here 


3  MISSION  PLANNING  ALGORITHMS 


In  this  section,  two  new  algorithms  are  presented  for  solving  the  two-dimensional 
mission  planning  problem.  The  projection  algorithm  begins  with  a  set  of  line 
segments  linking  M  to  T.  If  the  segments  intersect  any  of  the  threats,  then 
by  projecting  the  center  of  the  threat  circles  onto  the  nearest  line  segment, 
new  line  segments  can  be  constructed  which  usually  avoid  the  threats.  The 
circumscribed  triangle  algorithm  creates  paths  around  the  threats  by  creat¬ 
ing  triangles  which  enclose  the  threats.  The  paths  are  constructed  by  linking 
the  edges  of  these  special  trieingles. 

3.1  The  Projection  Algorithm 

The  projection  algorithm  begins  with  one  or  more  line  segments  linking  M 
(the  missile)  to  T  (the  target).  If  this  is  feasible,  then  the  algorithm  termi¬ 
nates.  Otherwise,  additional  break  points  are  created  to  guide  the  missile 
around  the  threats.  The  line  segment  [M,  T]  illustrated  in  Figure  2  has  been 
deflected  around  the  threat  by  adding  a  break  point  between  M  and  T . 

Let  Xi  denote  the  center  of  the  threat  circle  with  radius  ri .  Let  Z  be  the 
projection  of  Xi  onto  the  line  segment  [M,T].  Let  PF  be  a  new  point  ob- 


tained  by  moving  from  Xi  in  the  direction  Z  —  to  a  point  slightly  beyond 
the  circle.  The  new  path  is  [M,  W]  [W^,T]. 


Figure  2  about  here 


Deflecting  a  path  around  a  threat  can  result  in  the  new  path  passing 
through  a  different  threat  as  illustrated  in  Figure  3.  This  requires  another 
application  of  the  projection  procedure.  By  continuing  in  this  way,  we  Ccin 
in  most  cases  obtain  a  set  of  line  segments  that  link  M  to  T  and  avoid  all  of 
the  threats. 


Figure  3  about  here 


The  basic  algorithm  is  called  projector  &nd  the  algorithm  is  defined  as  follows; 


Algorithm :  projector 
Inputs  : 

M  -  missile  location 

T  -  target  location 

K  -  number  of  threats 

Xi,  ...,Xk  -  center  points  of  the  threats 

J*!)  *  ra-dii  of  the  threats 

L  -  number  of  segments  in  the  initial  path 

Ii, Jl+i  -  the  initial  path 


Outputs  : 

S  -  number  of  segments  in  solution 

Vi,  ...jYs+i  -  the  points  in  the  solution  path 


begin 

Ys  <—  Is  for  s  =  1, i  +  1; 
s  1,  fc  <-  1,  5  <-  i; 
while  (s  <  S)  do 
begin 

Z  <-  project 

if  {Z  is  on  the  line  segment  [y^Yi+i])  then 


begin 


if  (distance  (Z,  Xk)  <  rk)  then 
begin 

+  l.l(ri)(^  -  Xit)/distaiice  {Z,Xk) 
for  i  =  5  +  1  to  s  +  1  do  <— 

1^+1  <—  5  <—  5  +  l,s  <—  l,fc  <—  0; 

end; 
end; 

fc  <—  fc  +  1; 

if  [k  >K)  then  A:  <—  1,  s  s  +  1; 
end; 
end; 

procedure  distance  {U,V) 
begin 

return  the  Euclidean  distance  from  U  to  V; 
end; 
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procedure  project  {R,  U,  V) 
begin 

return  the  projection  of  R  onto  the  line  defined  by  points  U  and  V ; 
end; 


It  is  possible  for  Z  to  equal  Xk  so  that  Z  —  Xk  =  0.  When  this  occurs, 
we  perturb  Z  slightly  in  the  y-coordinate  direction.  Also,  since  the  threats 
overlap,  we  sometimes  obtain  a  path  with  zig-zagging  which  intersects  with 
one  or  more  circles.  At  the  conclusion  of  project,  we  call  a  scrubber  procedure 
which  removes  points  which  appear  inside  a  circle.  This  is  illustrated  in 
Figures  4  and  5. 


Figures  4  and  5  about  here 


To  obtain  a  path,  we  initialize  with  a  set  of  values,  rim  pro¬ 

jector  followed  by  the  scrubber.  We  compare  the  path  obtained  with  the 
shortest  feasible  path  found  so  far  and  retain  the  shortest  feasible  path.  The 


middle  or  direct  path  is  obtained  by  setting  7i  to  M,  I2  to  T,  and  X  to  1. 
The  high  path  attempts  to  find  a  path  such  that  all  threats  are  below  this 
path.  For  this  path,  I2  is  set  to  the  center  of  the  circle  first  encountered 
when  we  sweep  clockwise  a  line  segment  anchored  at  M  from  the  vertical  po¬ 
sition.  Let  Xk  denote  this  point  and  set  I2  to  Xh-  This  can  be  accomplished 
by  examining  the  slopes  of  the  line  segments  [M,  Xk\  and  selecting  the  one 
with  the  largest  slope.  The  point  J3  is  set  to  the  Xk  first  encountered  when 
we  sweep  counter-clockwise  a  line  segment  anchored  at  T  from  the  vertical 
position.  Finally,  74  is  set  to  T  and  L  is  set  to  3.  These  starting  segments 
generally  result  in  a  path  with  every  threat  below  the  path.  A  low  path  can 
be  constructed  in  a  similar  manner  by  modifying  the  selection  routine  for  I2 
and  I3.  A  middle  high  path  is  generated  by  selecting  I2  to  be  on  the  line 
segment  [M,Xh].  A  middle  low  path  is  constructed  in  a  similar  manner. 

3.2  The  Circumscribed  Triangle  Algorithm 

The  circumscribed  triangle  algorithm  starts  at  M  (the  missile  location)  and 
constructs  line  segments  by  moving  from  a  given  segment  point  tangentially 
to  the  circumference  of  nearby  threat  regions,  until  a  suitable  straight-line 
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path  from  a  segment  point  to  T  (the  target  location)  is  found.  This  basic 
strategy  is  embedded  in  a  branch- and-bound  framework,  so  that  a  minimal 
distance  path  of  this  restricted  type  is  obtained. 

In  this  algorithm  we  allow  the  missile  to  intersect  a  point  on  the  circum¬ 
ference  of  a  threat  region  but  do  not  allow  the  missile  to  fly  through  any 
interior  point  of  a  threat  region.  Further,  we  allow  a  limit  to  be  placed  on 
the  turn  angle  at  any  segment  point,  since  a  missile  is  unlikely  to  be  making  a 
sharp  turn  of,  say,  90®  or  more  at  such  a  point.  (Given  segments  [Vf-i,  Yi]  and 
[Yi,  Yi+i],  the  turn  angle  at  segment  point  Yi  is  defined  to  be  the  complement 
of  angle  y^-xYlYi+i.) 

Let  d{E,  F)  denote  the  Euclidean  distance  between  points  E  and  F.  The 
basic  step  in  proceeding  from  segment  point  Yi  tangentially  to  a  (nearby) 
threat  region  Rk  =  {s  ’•  d(Xk,s)  <  rjfe}  involves  constructing  the  circum¬ 
scribing  triangle  on  Rk  having  one  vertex  at  Yj  (see  Figure  6). 


Figure  6  about  here 


The  tangential  points  in  Figure  6  are  Ti ,  T2,  and  T3  and  the  vertex  points 


are  Yi,  Vi,  and  V2.  Sides  YiVi  and  ¥{¥2  are  the  equal  sides  of  the  (isosceles) 
circumscribing  triangle.  Tangential  point  T3  is  easily  computed  since  it  lies 
on  the  line  through  Yi  and  Xk  with  d{Xk,T3)  =  Vk-  To  compute  the  other 
points  we  make  use  of  a  result  from  elementary  geometry  which  states  that 
the  radius  r  of  a  circle  inscribed  in  a  triangle  with  sides  a,  6,  and  c  is  given 

by 

\/s{s  —  a)(s  —  b){s  —  c) 
r  =  — - 

s 

where  5  =  |(a  +  6  +  c). 

Letting  a  =  d{Yi,Vi),  b  =  d{Yi,V2),  and  c  =  d(Vi,V2)  with  a  =  6,  a 
one- dimensional  search  on  the  value  of  c  can  be  used  to  produce  a  value  of 
r  =  Tfc.  From  this  c  value,  points  Vi  and  V2  can  be  determined.  Points  Ti 
and  T2  can  then  be  computed  as  the  nearest  points  to  Xk  on  YiVi  and 
respectively. 

We  intend  to  mahe  use  of  these  tangential  lines  to  construct  the  next 
segment,  if  possible.  Before  using  points  on  either  tangential  line,  we  first 
make  sure  that  the  respective  turn  angles  at  Yi  would  be  within  the  turn 
angle  limit. 

The  points  Vi  and  V2  have  a  special  property.  Each  is  the  first  point  on  the 


tangential  line  from  Yi  through  itself  which  has  a  complete  180°  field  of  view 
(ignoring  turn  angle  constraints)  of  the  side  of  Rk  opposite  Yi.  By  contrast, 
tangent  points  Ti  and  Tj  each  have  less  than  a  90°  field  of  view  (ignoring 
turn  angle  constraints)  of  that  side  (see  Figure  7).  With  this  in  mind  Vi  (or 
V2)  would  seem  to  be  a  good  choice  for  the  end  point  of  a  next  segment  in 
our  path  construction.  We  actually  attempt  to  use  a  point  between  Ti  and 
Vi  (or  between  Tj  and  V2)  for  the  next  end  point.  In  general,  such  a  point 
can  be  represented  as  a  convex  combination  Vi+i  =  (1  —  oc)Ti  +  aVi  (or 
=  (1  -  a)T2  +  ccV2  ),  for  some  0  <  a  <  1.  In  order  for  such  a  point 
ii+i  to  be  usable,  we  must  insure  that  the  segment  is  clear.  (A  line 

segment  is  defined  to  be  clear  if  it  does  not  intersect  the  interior  of  any  of 
the  threat  regions.) 

We  first  search  for  a  minimum  value  of  a  with  0  <  a  <  1  such  that 
[y;-,  (1  - a)Ti  +  aVi]  (or  [Vi,  (1  - a)T2  +  aV2] )  is  clear  and  [(1  - oc)Ti  +  aVuT] 
(or  [(1  —  oc)T2  +  aV2,T]  )  is  clear  without  exceeding  the  turn  angle  limit  at 
(1  -  a)Ti  +  aVi  (or  (1  -  a)T2  +  aV2  ).  If  such  a  point  is  found  we  have  a 
feasible  completion  by  adding  the  two  segments  [Fijlf+i  =  (1  —  oc)Ti  +  qVi] 
and  [y;+i  =  (1  -  a)Tk  +  aVuYi+2  =  T]  (or  [Y^Y+i  =  (1  -  a)T2  +  aV2]  and 
[Yi+r  =  {l-a)T2  +  aV2,Yi+2  =  T]). 
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If  such  a  point  leading  to  a  feasible  completion  cannot  be  found,  it  would 
seem  natural  to  use  Vi  (or  V2)  by  default,  due  to  its  wide  field  of  view. 
However,  in  practice  this  leads  to  longer  paths  and  we  find  it  better  to  use 
a  point  between  T\  and  Vi  (or  between  T2  and  Vi)  for  the  next  end  point. 
We  have  parameterized  this  by  defining  a  backoff  percentage  0  which  specifies 
that,  in  this  case,  Yi  is  to  be  chosen  as  the  convex  combination  (1— ;0)ri+^Vi 
(or  (1  —  /5)T2  +  /3V2  )•  Experimentally,  we  have  found  ^  |  to  be  a  good 

choice. 


Figure  7  about  here 


The  branch-and-bound  algorithm  we  have  developed  based  on  this  cir¬ 
cumscribed  triangle  approach  is  given  as  follows: 

Algorithm  :  circumscribe 
Inputs  : 

M  -  missile  location 

T  -  target  location 


K 


-  number  of  threats 


Xi,...,Xk 

ri,...,rK 

Z 


Outputs  : 


D 

S 


Li, Ls 
Constructs  : 

B 


Q 

d{E,F) 

L(E,F) 

R(C,r) 

S{C,r) 

I 


Yi 

Si 


-  center  points  of  the  threats 

-  radii  of  the  threats 

-  turn  angle  limit  (degrees) 

-  fall  back  parameter  (0  <  ^  <  1) 

-  length  of  best  constructed  solution  path 

-  number  of  segments  in  best  constructed  solution  path 

-  segments  of  best  constructed  solution  path 

-  pointer  to  last  point  in  best  constructed  solution  path 

-  branch-and-bound  queue 

-  distance  from  E  to  F 

-  line  segment  from  E  to  F 

-  threat  region  with  center  C  and  radius  r 

-  interior  of  threat  region  with  center  C  and  radius  r 

-  number  of  points  constructed 

-  constructed  point 

-  pointer  to  point  previous  to  Yi 


Di 

Si 

Ei 


Tx,T2,T2 


Vi,y2 


-  path  distance  prior  to  Yi 

-  path  segments  prior  to  Yi 

-  set  of  threat  points  eligible  for  moving  away  from  point  Yi 

-  constructed  tangent  points  of  circumscribled  triangle 

-  constructed  vertices  of  circumscribled  triangle 


begin 

jD  <—  oo;  S  *-0]  B  <—  0; 
if  L{M,T)r\S{Xi,ri)  =  ^  for  i  =  1..A' then 
D  ^  d{M,Ty,  5^1;  h  ^  [M,T]; 
else 
begin 


/  <— 1;  Yi  <— M;  Bi  <— 0;  <— 0;  5i  <— 0;  Ei  *— {Ri,...,Rk}‘, 


put  {Yx,Bi,Di,Si,Ei)  on  Q] 


while  (5  7^  $  do 
begin 


remove  some  (Y^-,  A,  Si,  Ei)  from  Q; 


A; 
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select  some  Rj  from  jB,-;  E{  *—  E{  \  Rf, 

if  Ei  ^  $  then  put  {Yi,Bi,Di,Si,Ei)  back  on  Q] 

form  the  triangle  circumscribing  Rj  having 

its  tangent  points  Ti  G  YiVi,T2  €  YV2,Tz  G  ViV2‘, 

for  fc  =  1  to  2  do 

begin 

Um-YbYiTk<Z  and 

L{Yi,Tk)r\S{Xi,ri)  =  ^  for  l=l..K  and 
Di  +  d{Yij  Tk)  +  d{Tk,  T)  <  D  then 
begin 

search  for  a  point  U  on  L{Tk,  Vk)  nearest  to  Tk  satisfying: 
m-YUT<Z  and 

L{Yi,U)nS{Xi,ri)  =  il>  for  I  =  1..K  and 
L{UJ)r\S{Xuri)  =  ^  for  l  =  l..K] 
if  such  a  point  U  exists  then 
begin 

if  A  +  d{Yi,  U)  +  d(  A  T)<D  then 
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begin 

/  ^7  +  1; 

Yi^U\  Bi^i;  Di^Di  +  d{Yi,Uy,  5/ ^5^  +  1; 

m  i —  7j  I  < —  7-1-1’, 

Yj  *— T',  Bj  ^  m\  Dj  <—  Djn  -{-  d{U,T)\  Si  *—  Si  2; 

7)^7);;  S^  Si\  B^I] 
end; 
else 
begin 

if  Ei  ^  $  then 
begin 

7^  ^  (1  -  -h 

if  L{Yi,U)r\S{Xi,ri)  =  ^  for  /=1..7i:then 
begin 
7<-7-l-l; 

Yi^U-,  Bi^y  Di^Di^d{Yi,Uy  Si^Si  +  V,  Ei  ^  Ei 
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put  {Yi,Bi,Di,Si,Ei)  on  Q; 
end; 
end; 


end; 

end; 

end; 

end;  • 

end; 

i  ir-  S]  k  *—  B] 

while  fc  0  do  j  *— Bk\  ?•*— i  — 1;  k  *— j; 

end; 


Figure  8  shows  an  optimal  path  through  a  region  containing  three  threats 
generated  by  the  circumscribed  triangle  algorithm.  Figure  9  shows  additional  • 

partial  paths  produced  by  the  branch-and-bound  process  contained  in  the 
algorithm.  The  actual  implementation,  whose  results  are  reported  in  Section 
IV,  actually  malces  two  applications  of  the  algorithm,  the  second  of  which  has 
the  roles  of  missile  and  taxget  reversed.  Of  course,  any  path  obtained  with 
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missile  and  target  reversed  would  itself  need  to  be  reversed.  The  shortest  of 
the  two  paths  generated  is  selected  eis  the  best  path. 


Figures  8  and  9  about  here 


4  EMPIRICAL  ANALYSIS 

Using  a  grid  size  of  780  x  700,  we  randomly  generated  five  groups  of  test 
problems.  Each  group  had  five  SAM  sites  with  radii  randomly  generated  on 
the  interval  [50,  150].  The  four  problems  in  each  group  had  identical  SAM 
sites,  but  different  missile  and  target  locations  with  the  missiles  located  on 
the  left  boundary  and  targets  located  on  the  right  boundary.  The  missile 
and  target  locations  were  randomly  selected  and  were  different  for  each  of 
the  twenty  test  problems.  The  input  for  each  problem  consists  of  seven  points 
(the  missile  location,  the  target  location,  five  SAM  sites)  and  five  radii. 

A  graphical  system  was  developed  using  Tcl  and  Tk  (see  Welch  [6]).  Tcl 
stands  for  Tool  Command  Language  and  Tk  is  a  toolkit  for  window  pro¬ 
gramming.  The  system  allowed  us  to  display  the  problem  graphically  and 


construct  a  path  manually  by  simply  pointing  and  clicking.  This  system  was 
used  to  manually  construct  solutions  for  each  of  the  twenty  test  problems. 
Both  the  projection  and  the  circumscribed  triangle  algorithms  were  also  ap¬ 
plied  to  each  of  the  test  problems.  The  results  are  summarized  in  Table  1. 
The  manual  paths  tended  to  steer  clear  of  the  SAM  sites  and  use  only  a 
few  segments  at  the  penalty  of  substantially  longer  paths.  The  projection 
algorithm  generated  relatively  short  paths  at  the  expense  of  substantially 
more  segments.  The  circumscribed  triangle  algorithm  worked  best,  using  a 
small  number  of  segments  and  producing  very  short  paths.  The  paths  ob¬ 
tained  manually  and  by  the  projection  algorithm  may  be  found  in  Figures  10 
through  14  in  Appendix  A.  The  paths  obtained  by  the  circumscribed  triangle 
algorithm  may  be  found  in  Figures  15  through  19  in  Appendix  B. 


Table  1.  Comparison  of  Three  Procedures  for  the  Mission  Planning  Problem 


Manual 

Projection  Alg 

Circumscribed  Triangle  Alg 

Distance 

Segments 

Distance 

Segments 

Distance 

Segments 

2 

959 

3 

957 

3 

3 

853 

5 

853 

3 

998 

4 

948 

5 

853 

4 

3 

946 

4 

945 

3 

826 

2 

798 

2 

795 

2 

934 

4 

881 

4 

878 

4 

7 

2 

1043 

4 

1042 

3 

8 

813 

2 

807 

2 

805 

2 

9 

2 

784 

2 

782 

2 

10 

3 

1054 

4 

1049 

3 

11 

3 

852 

3 

847 

3 

12 

839 

2 

823 

2 

820 

2 

13 

822 

2 

811 

2 

810  ! 

2 

14 

■ukliH 

3 

946 

5 

945 

3 

15 

3 

961 

5 

966 

3 

16 

877 

3 

839  . 

3 

837 

3 

17 

911 

3 

864 

4 

865 

3 

18 

3 

1014 

4 

1009 

3 

19 

4 

974 

4 

900 

4 

20 

3 

855 

3 

862 

2 

Totals 

18,984 

56 

18,012 

70 

17,820 

57 

1.07 

0.98 

1.01 

1.23 

1.0 

1.0 

5  SUMMARY  AND  CONCLUSIONS 


This  manuscript  presents  two  algorithms  for  the  mathematical  problem  of 
finding  a  path  consisting  of  line  segments  from  a  given  origin  to  a  given 
destination  which  does  not  pass  through  any  of  a  set  of  threats  represented  by 
circles.  Neither  algorithm  requires  the  use  of  a  grid  imposed  on  the  military 
theatre  which  is  required  by  the  competing  algorithms  of  which  we  are  aware. 

The  projection  algorithm  initially  ignores  the  threats  and  constructs  a 
path  from  the  origin  to  the  destination.  Let  [Yj,  Yj+i]  denote  the  line  seg¬ 
ment  defined  by  endpoints  Yj  and  Y^+i.  If  [Yj,  Yj+i]  passes  through  a  threat 
whose  center  is  X,  then  X  is  projected  onto  [Yi,Yi+i]  to  obtain  Z.  It  next 
finds  a  point  W  outside  the  threat  on  the  line  defined  by  Z  and  X.  The 
original  segment  [YijYj+i]  is  replaced  by  two  segments  [YiiVY"]  and  [WjYj+i]. 
Continuing  in  this  way  usually  leads  to  a  feasible  path. 

The  circumscribed  triangle  algorithm  constructs  triangles  around  the 
threats  and  selects  break  points  on  either  the  upper  or  lower  edge  of  these 
circumscribed  triangles.  Hence,  from  the  origin  there  are  two  potential  paths 
around  a  threat  (the  upper  path  and  the  lower  path).  From  each  of  these 
paths  there  may  be  additional  threats  from  which  we  consider  both  the  up- 


per  path  and  the  lower  path.  This  combinatorial  feature  of  the  strategy  is 
handled  in  a  branch-and-bound  framework. 

In  an  empirical  analysis  with  twenty  randomly  generated  test  problems, 
we  found  that  the  circumscribed  triangle  algorithm  consistently  produced 
paths  which  were  superior  to  those  obtained  either  manually  or  by  the  pro¬ 
jection  algorithm.  Also,  our  software  implementation  is  so  fast  that  this 
algorithm  can  be  run  on  a  PC  in  a  real-time  environment. 
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APPENDIX  A 


SOLUTIONS  USING  PROJECTION  ALGORITHM 

This  appendix  presents  a  graphical  display  of  the  twenty  randomly  generated 
test  problems  along  with  the  mission  obtained  manually  and  the  mission 
obtained  by  the  projection  algorithm.  The  black  line  is  the  mission  obtained 
using  projector  and  the  gray  line  is  the  mission  obtained  manually. 


Figures  10  to  14  about  here 


APPENDIX  B 


SOLUTIONS  USING  THE  CIRCUMSCRIBED 
TRIANGLE  ALGORITHM 


This  appendix  presents  a  graphical  display  of  the  solutions  produced  by  the 
circumscribed  triangle  algorithm. 


Figures  15  to  19  about  here 


Figure  1;  Illustration  of  a  Feasible 
Five  SAM  Sites 
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Figure  9;  Partial  Paths  Considered  by  the  Circumscribed  Triangle  Algorithm 
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Figure  10:  Display  of 
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1,  2,  3,  and  4 


Figure  11:  Display  of  Problems  5,  6,  7,  and  8 
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Problem  9 


Problem  10 


Figure  12:  Display  of  Problems  9,  10,  11,  and  12 
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Figure  16:  Display  of  Problems  5,  6,  7,  and  8 
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Figure  17:  Display  of  Problems  9, 10,  11,  and  12 
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Figure  19;  Display  of  Problems  17, 18,  19,  and  20 
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Abstract 


By  representing  threats  as  circles,  the  problem  of  mission  planning  for  a  cruise 
missile  can  be  viewed  as  a  computational  geometry  problem  in  which  we  seek 
a  path  composed  of  line  segments  from  a  launch  site  to  a  target  site  which 
skirts  selected  threats.  In  a  previous  investigation,  we  presented  an  effective 
algorithm  of  this  type,  called  the  circumscribed  triangle  algorithm.  In  this 
investigation,  we  extend  the  circumscribed  triangle  algorithm  to  incorporate 
a  side  constraint  on  the  probability  of  a  successful  traversal.  This  allows  the 
missile  to  take  short  cuts  through  threat  areas,  as  long  as  the  probability  of 
success  meets  some  user-specified  criteria.  In  an  empirical  study  on  twenty 
randomly  generated  test  problems  using  probabilities  of  90%,  75%,  and  60%, 
we  found  that  the  new  algorithm  works  extremely  well. 
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1  INTRODUCTION 


The  cruise  missile  has  recently  become  one  of  the  Navy’s  most  important 
weapon  systems.  During  recent  conflicts  in  the  Persian  Gulf,  ship  and  subma¬ 
rine  launched  cruise  missiles  have  been  used  extensively  with  great  success. 
Targets  which  may  be  hundreds  of  miles  inland  can  be  attacked  by  these 
Navy  weapons  stationed  in  relatively  safe  positions  offshore.  In  addition, 
cruise  missiles  can  be  used  to  disable  SAM  sites,  making  manned  aircraft 
strikes  substantially  safer. 

A  mission  for  a  cruise  missile  is  defined  as  a  path  from  a  launch  area  to 
the  target.  The  job  of  the  mission  planner  is  to  determine  the  path  that  the 
missile  will  follow  and  program  the  missile  with  the  appropriate  waypoints 
and  any  TERCOM  maps  that  will  be  used.  A  mission  does  not  have  to 
include  TERCOM  maps,  but  they  can  be  used  as  part  of  the  navigational 
package.  Some  missions  involve  only  GPS  waypoints. 

Using  current  technology,  missions  are  developed  using  a  two-step  process. 
In  the  first  step,  a  mission  planner  uses  a  two-dimensional  map  to  manually 
select  a  path  from  the  launch  site  to  the  target.  A  software  system  enhances 
the  path  to  include  the  vertical  dimension  and  develops  an  estimate  of  the 


probability  that  the  mission  will  be  successful.  This  process  may  be  repeated 
several  times  until  the  mission  planner  is  satisfied  with  the  mission  plan. 
Since  a  cruise  missile  costs  approximately  $650,000,  the  probability  of  success 
of  a  given  mission  is  an  important  consideration  for  Navy  management. 

Several  research  groups  have  attempted  to  develop  automatic  methods  to 
replace  the  manual  part  of  the  two-step  approach.  Most  of  these  groups  begin 
by  placing  a  grid  over  the  combat  theatre  and  then  applying  a  modification 
of  Dijstra’s  [2]  algorithm  to  obtain  a  path  from  the  origin  to  the  destination 
(see  [1,  3,  5,  6,  8,  9]).  Unfortunately,  the  grid  can  become  very  large  and 
consequently  the  software  can  be  fairly  slow.  In  addition,  the  resulting  path 
can  involve  numerous  line  segments. 

In  [4],  the  authors  developed  an  effective  algorithm  for  mission  planning 
that  does  not  involve  the  use  of  a  grid.  The  threats  are  modelled  as  circles 
on  a  two-dimensional  map,  and  the  objective  is  to  obtain  a  shortest  path 
from  a  given  origin  to  a  given  destination  which  does  not  pass  through  any 
threat  circle.  The  circumscribed  triangle  algorithm  constructs  line  segments 
by  moving  from  a  given  segment  point  tangentially  to  the  circumference  of 
nearby  threat  regions,  until  a  suitable  path  from  origin  to  destination  is 
found.  In  an  empirical  study,  the  algorithm  always  obtained  a  very  good 
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mission  plan. 

The  objective  of  this  investigation  is  to  incorporate  a  probability  con¬ 
straint  into  the  circumscribed  triangle  algorithm.  This  allows  for  more  re¬ 
alistic  mission  plans  which  fly  through  threat  regions  on  their  path  to  the 
target  while  maintaining  a  user-specified  probability  of  success.  The  benefit 
of  this  extension  is  determined  empirically  on  a  test  suite  of  twenty  randomly- 
generated  problems. 

2  THE  MISSION  PLANNING  PROBLEM 

In  this  work,  a  threat  is  represented  by  a  circle  with  a  given  center  X  and  a 
given  radius  r.  Let  Z  denote  the  point  along  the  mission  plan  at  which  the 
missile  is  nearest  X.  If||Z  —  X||  >r,  then  the  threat  at  X  is  too  far  away  to 
harm  the  missile  and  the  probability  of  a  successful  traversal  is  defined  to  be 
1.  Otherwise,  some  function  p{Z,X,r)  is  used  to  determine  the  probability 
of  a  successful  traversal.  When  there  are  multiple  threats,  the  probability  of 
a  successful  traversal  is  the  product  of  the  probabilities  along  each  segment, 
taking  into  account  all  of  the  threats. 

The  problem  of  finding  a  feasible  mission  plan  which  incorporates  a  prob- 


ability  constraint,  can  be  defined  mathematically  as  follows: 


Given  a  target  probability  P,  two  points  M  and  T  €  B?,  and  K 
circles,  with  centers  Xi,...,Xk  €  R'^  and  radii  rj, ....tk, 
find  a  set  of  line  segments  [YuY^],  [Y^^Yi],...,  [Is,i5+i], 
with  Yi  =  M  and  =  P,  such  that  the  probability  of  a 
successful  traversal  is  at  lecist  P. 

The  quality  of  a  mission  plan  is  measured  by  the  length  of  the  path  and  the 
number  of  path  segments.  The  best  path  is  the  straight  line  from  M  to  T.  If 
this  is  not  feasible  (i.e.  the  probability  of  a  succcessful  traversal  is  less  than 
P),  then  we  seek  a  short  path  having  only  a  few  segments. 

3  THE  ALGORITHM 

The  probabilistic  incursion  triangle  algorithm  allows  the  missile  to  fly  through 
the  interior  of  a  threat  region  and  employs  a  probability  model  to  compute 
the  probability  of  survival  for  a  path  from  the  missile  location  M  to  the  tar¬ 
get  location  T  which  consists  of  line  segments  which  may  either  be  tangent 
to  a  circular  threat  region  or  cut  through  the  region  (an  incursion). 

The  algorithm  operates  within  the  same  branch-and-bound  algorithmic 
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structure  used  in  the  circumscribed  triangle  algorithm  (see  Section  3.2  of 
[4]).  The  major  differences  are  (1)  when  attempting  to  construct  a  new 
segment  by  approaching  a  threat  region  from  a  previous  segment  endpoint, 
several  equilateral  triangles  passing  through  the  threat  region  in  addition  to 
the  circumscribed  triangle  may  be  considered,  and  for  each  such  triangle, 
potentially  two  points  are  added  to  the  branch-and-bound  queue,  and  (2)  in 
searching  for  a  final  segment  linking  a  previous  segment  endpoint  to  the  tar¬ 
get,  that  segment  may  intersect  threat  regions  (with  a  consequent  reduction 
in  survival  probability). 

3.1  The  Probability  Model 

For  threat  region  k  with  center  Xk  and  radius  rk,  the  probability  model 
assumes  that  the  probability  of  survival  of  a  line  segment  approaching  a 
threat  region  is  a  function  only  of  the  minimum  distance  from  the  segment 
to  Xk-  Further,  the  model  assumes  independence  so  that  the  overall  survival 
probability  for  a  path  is  the  product  of  the  survival  probabilities  for  each 
segment  in  the  path. 

For  developmental  purposes  we  have  hypothesized  a  crude  probability 
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function  which  is  a  quadratic  of  the  form 


Pr{d)  =  ad^  +  c  , 

where  for  threat  region  k, 

Pr{rk)  =  .99  and  Pr(0)  =  .50  . 

3.2  Incursion  Triangles 

Figure  6  illustrates  a  circumscribed  triangle  YiViV2  and  an  incursion  trian¬ 
gle  YiQiQ2  for  a  common  threat  region  Rk-  In  the  circumscribed  triangle 
algorithm,  with  respect  to  the  segments  YiVi  and  YiV2,  the  line  segments 
TiVi  and  T2V2  are  searched  first  for  a  clear  completion  segment  to  the  tar¬ 
get,  before  the  points  Vi  and  V2  are  added  to  the  branch-and-bound  queue. 
This  also  occurs  in  the  probabilistic  incursion  triangle  algorithm,  but  the 
completions  can  intersect  threat  regions  if  the  total  survival  probability  is  at 
least  P  and  the  segments  YiVi  and  5^1^  will  each  have  survival  probability 
.99.  Also,  in  the  probabilistic  incursion  triangle  algorithm,  with  respect  to 
the  segments  KQi  and  YiQ2,  the  line  segments  PiQi  and  P2Q2  are  searched 
first  for  a  completion  to  the  target  with  total  survival  probability  at  least  P, 
before  the  points  V\  and  V2  are  added  to  the  branch-and-bound  queue. 


The  number  of  incursion  triangles  considered  in  approaching  a  threat 
region  from  a  previous  segment  endpoint  depends  on  the  accumulated  number 
of  segments  and  an  adjustable  parameter  S  which  is  roughly  interpreted  as 
the  expected  number  of  segments  in  an  optimal  path.  In  the  testing  described 

A 

in  Section  4,  S  was  set  at  3.  If  Si  is  the  number  of  segments  in  a  partial 
path  from  the  missile  location  to  segment  endpoint  Yi  and  P,-  is  the  survival 
probability  of  this  partial  path,  incursion  triangles  with  individual  segment 
survival  probabilities  of 

.99{P/Pi)  ,  •99(  ^  ■  •99(  ^  ) 

will  be  considered,  unless  these  probabilities  are  very  close  to  .99. 

4  EMPIRICAL  ANALYSIS 

The  twenty  test  problems  presented  in  [4]  were  solved  using  P  =  90%, 
P  =  75%,  and  P  =  60%,  and  the  results  axe  summarized  in  Table  1.  The 
deterministic  algorithm  is  the  circumscribed  triangle  algorithm  described  in 
[4].  Note  that  setting  P  —  60%  subsantially  reduces  the  number  of  way- 
points,  but  only  reduces  the  total  distance  by  4%.  The  paths  obtained  by 
the  deterministic  algorithm  and  the  new  algorithm  with  P  =  60%  may  be 


found  in  Appendix  A.  For  most  problems,  the  two  paths  are  similar,  with 
the  new  algorithm  taking  short  cuts  through  one  or  more  threats.  Problem  7 
illustrates  a  case  where  the  new  algorithm  selects  a  completely  different  path. 
For  this  problem,  the  deterministic  algorithm  produces  a  path  in  which  all 
threats  are  above  the  path,  while  the  new  algorithm  uses  a  mixed  strategy. 
The  deterministic  strategy  is  unable  to  use  a  single  straight  line  for  any  of 
these  test  problems,  whereas  the  probabilistic  algorithm  finds  this  to  be  best 
for  several  of  them.  Based  on  these  results,  we  believe  that  this  approach 
could  be  used  in  the  development  of  a  grid-free  auto-router,  incorporating  a 
probabilistic  side  constraint. 


Table  1.  Comparison  of  Probabilities  for  the  Mission  Planning  Problem 


Prob 

Deterministic 

Probabilistic  Algorithm 

Algorithm 

90% 

75% 

60% 

Dist 

Segmts 

■BBl 

Segmts 

Dist 

Segmts 

Dist 

Segmts 

1 

957 

3 

955 

3 

952 

3 

949 

3 

2 

853 

3 

849 

3 

842 

3 

834 

3 

3 

853 

4 

851 

3 

843 

3 

803 

3 

4 

945 

3 

939 

3 

925 

3 

905 

3 

5 

795 

2 

794 

2 

789 

2 

783 

2 

6 

878 

4 

876 

4 

868 

4 

857 

4 

7 

1042 

3 

1040 

3 

1031 

3 

1012 

3 

8 

805 

2 

805 

2 

804 

1 

804 

1 

9 

782 

2 

782 

2 

780 

2 

780 

1 

10 

1049 

3 

1044 

3 

1031 

3 

1012 

3 

11 

847 

3 

845 

3 

838 

3 

835 

2 

12 

820 

2 

819 

2 

817 

2 

815 

1 

13 

810 

2 

809 

2 

807 

2 

807 

1 

14 

945 

3 

942 

3 

861 

3 

809 

2 

15 

966 

3 

962 

3 

950 

3 

780 

2 

16 

837 

3 

837 

3 

835 

3 

835 

1 

17 

865 

3 

857 

3 

848 

3 

825 

3 

18 

1009 

3 

1005 

3 

994 

3 

977 

3 

■a 

900 

4 

895 

3 

882 

3 

857 

3 

862 

2 

851 

3 

837 

2 

820 

3 

Totals 

17820 

57 

kb 

56 

17534 

54 

17099 

47 

Scaled 

1.0 

|■QS| 

0.982 

0.984 

0.947 

0.960 

0.825 
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APPENDIX  A 


SOLUTIONS 

This  appendix  presents  a  graphical  display  of  the  twenty  test  problems  along 
with  the  solutions  obtained  by  the  deterministic  algorithm  and  the  prob¬ 
abilistic  algorithm.  The  black  line  is  the  solution  obtained  with  a  P  set  to 
60%  and  the  gray  line  is  the  mission  obtained  by  the  deterministic  algorithm. 
The  plots  were  obtained  using  Tcl  and  Tk  Solver  [7],  which  is  part  of  our 
experimental  computational  package. 
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Figure  2:  Display  of  Problems  5,  6,  7,  and  8 
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Figure  3:  Display  of 
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Problem  10 


9,  10,  11,  and  12 
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Figure  5:  Display  of  Problems  17,  18,  19,  and  20 
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Figure  6:  Circumscribing  Triangle  and  Incursion  Triangle  with  Vertex  Yi  for 
Threat  Region  Rk 
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